Vortex states and Majorana fermions in spin-orbit coupled semiconductor-superconductor hybrid 

structures 
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We study the energy spectrum of a vortex core in a two-dimensional semiconductor with Rashba spin-orbit 
interaction and proximity-coupled to a conventional superconductor and a ferromagnetic insulator. We perform 
self-consistent calculations using the microscopic tight-binding Bogoliubov-de Gennes method on a lattice and 
confirm the existence of Majorana fermions in the non-trivial topological phase. We also find two different 
topologically trivial superconducting phases, only differing in the type of vortex core structure they support 
and separated by a normal fermionic zero-energy excitation. Furthermore, we find an asymmetry in the energy 
spectrum with respect to both Zeeman splitting and vortex rotation direction and explain its physical origin. 

PACS numbers: 74.90.+n, 03.65.Vf, 03.67.Lx, 74.45.+C 



I. INTRODUCTION 

In condensed matter physics the topic of topology has a rich 
history. For several decades topology has played an impor- 
tant role in the study of phenomena such as nematics, do- 
main walls, and vortices, through the winding numbers of 
the corresponding order parameters.' ^ With the discovery of 
the quantum Hall effect^^ ^ a second type of topological the- 
ory emerged, in which topology in reciprocal rather than real 
space is important.^ Inspired by the topological concepts in- 
troduced through the quantum Hall effect, the existence of a 
new class of materials called topological insulators was re- 
cently theoretically predicted*^ ^ and subsequently experimen- 
tally verified.Q^ This has opened up a large research field, 
which also includes the very closely related concept of topo- 
logical superconductors." '^ 

On an initially unrelated note, Ettore Majorana introduced 
in 1937 a slightly modified version of the Dirac equation, 
where the particles are their own anti-particles. - The neutrino 
has been proposed to be a Majorana fermion, but experimen- 
tal investigations has so far failed to confirm the existence of 
any fundamental fermions of the Majorana type.''* More re- 
cently, Majorana fermions have been proposed to exist as ef- 
fective quasi-particle excitations on edges or in vortex cores 
of either spinless p + ip-wave superconductors or fractional 
quantum Hall systems with filling fraction v = \ .G^Qll From 
a theoretical point of view this provides an interesting play- 
ground, where concepts from both real and reciprocal topo- 
logical theories, in the form of vortices and quantum Hall-like 
physics, respectively, meets high-energy ph ysics in the long 
standing quest to discover Majorana fermions. E3EI1 In addition, 
and application-wise more relevant, the non-Abelian nature of 
the Majorana fermions makes them compelling candidates for 
implementing fault-tolerant topological quantum computation 
through braiding schemes.ESl 

Recently, an effective two-dimensional (2D) spinless p+ip- 
wave superconducting state has been proposed in both topo- 
logical insulators^** and ferromagnetic Rashba spin-orbit cou- 
pled semiconductors'^^' ^^ in proximity to a conventional s- 
wave superconductor. The maturit y of semiconductor tech- 
nology, large spin-orbit couplingp^'^^' and experimentally 



demonstrated superconducting proximity effecP^l rnake es- 
pecially the hybrid semiconductor-superconductor alternative 
very promising and possible experimental signatures of Majo- 
rana physics has already been reported in such nanowires.^" SI 
Although proposals exists for braiding operations in a network 
of ID nanowires,^^ vortices offer a more direct experimental 
route to braiding, important for both testing the non-Abelian 
nature of Majorana fermions as well as implementing compu- 
tational operations. 

While several studies based on a continuum model have 
established the existence of a Majorana mode in the voitex 
core in feiTomagnetic hybrid semiconductor-superconductor 
structures,'22E334| j^^ detailed microscopic self-consistent 
study exists. In this article we therefore study a micro- 
scopic model of a 2D Rashba spin-orbit coupled semicon- 
ductor with proximity-induced s-wave superconductivity and 
Zeeman exchange splitting. We perform self-consistent cal- 
culations using the microscopic tight-binding Bogoliubov-de 
Gennes method^^ on a lattice and obtain Majorana fermion 
solutions when the parameters are tuned such that the system 
enters the topologically non-trivial phase. In addition to the 
topologically non-trivial phase we find two topologically triv- 
ial superconducting phases, distinguished by their vortex core 
magnetization. The two phases are separated by a normal 
zero-energy vortex core mode, explicitly demonstrating that 
zero-energy states in these hybrid structures are not necessar- 
ily Majorana fermions. Moreover, we find a distinct asymme- 
try in the vortex core energy spectrum as function of either 
Zeeman splitting or vortex rotation direction. We show that 
this asymmetry is also present in the continuum model and 
provide a physical understanding of its origin. These results 
establish a remarkably rich behavior of the vortex core energy 
spectrum, including different zero-energy states. With the en- 
ergy spectrum being the primary experimental screening tool 
for identifying Majorana fermions, our results should provide 
a guideline in the search for vortex Majorana states. 



II. METHOD 

We consider a 2D semiconductor on a square lattice with 
Rashba spin-orbit interaction, spin-singlet s-wave supercon- 



ductivity and Zeeman splitting. This system can be realized 
in a thin Rashba spin-orbit coupled semiconductor layer sand- 
wiched between a conventional superconductor and a ferro- 
magnetic insulator, which by proximity-effect induce the s- 
wave superconductivity and Zeeman splittin g, see Fi g.fT] The 
Hamiltonian for this system can be written aJ2ll22|36| 
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(1) 




Here i and j are site indices on the square lattice, a is the 
spin index, and cl^{cic) is the electronic creation (annihi- 
lation) operator, t, jjl, V^, a and A are the nearest neigh- 
bor hopping, chemical potential, Zeeman splitting, Rashba 
spin-orbit interaction, and superconducting order parameter, 
respectively. For concreteness, we consider a lightly hole- 
doped semiconductor by setting /x = 4 and measure all en- 
ergies relative to the kinetic term, which we set to t == 1. We 
are interested in the properties of superconductor vortices in 
this system, for which the order parameter can be written as 
A(r, 6') = |A(r', 6')|e™^ around some vortex center r = 0, 
where n — ±.\ determines the vortex rotation direction. 

We solve Eq. ([TJ self-consistently within the Bogoliubov- 
de Gennes formulation.^S Due to the spin-orbit coupling and 
Zeeman field, a doubling of the number of degrees of freedom 
is necessary, resulting in two eigenstates for every ordinary 
electronic degree of freedom. The self-consistent calculation 
is carried out by specifying a vortex-like initial configuration 
^(0) _ e'"^. The superconducting pairing potential V^c is 
then introduced and used in subsequent steps to calculate 



E„<Q 



(2) 



where Ui,v<^(yv\i^ is the electron up (hole down) component on 
site i in the z/'th eigenstate and m is the iteration step. A('"+^) 
is then fed back into the calculation until a well converged A 
is obtained. It is thus the pairing potential Vgc together with 
vortex rotation direction n that self-consistently determines 
A. If not otherwise stated, calculations were done on a lattice 
of size 39 x 39. 

The Hamiltonian in Eq. (fill has several topologically non- 
trivial phases. For a lightly hole-doped semiconductor as con- 
sidered here, the condition for being inside the only experi- 
mentally accessible topologically non-trivial phase is^ 



FIG. 1: (Color online) Schematic setup with (a) thin semiconduc- 
tor layer with Rashba spin-orbit interaction into which .s-wave super- 
conductivity and Zeeman splitting are introduced by proximity effect 
from (b) a conventional superconductor and (c) a ferromagnetic in- 
sulator. We study (d) a superconducting vortex in this system. 



However, the superconducting order parameter is determined 
self-consistently and it will respond to changes in V^ and ap- 
proach zero for moderately to large V^ ■ At this point the super- 
conducting gap is destroyed and the topological band theory 
breaks down. As this will happen long before the upper bound 
in Eq. p) is violated, the condition for the non-trivial phase 
can be written as 



< lAI < |T4 



(4) 



|A| < |T4| < Va^2 + |a|2. 



(3) 



Having fixed t = 1, /i = 4, and replaced A with Vsc and 
n as external parameters there are four remaining parameters 
in the Hamiltonian in Eq. (fT). These are the strength of the 
Zeeman splitting (Vj), Rashba spin-orbit interaction (a), su- 
perconducting pairing potential (Kic) and vortex rotation di- 
rection (n). 



III. RESULTS AND DISCUSSION 

A. No Rashba Interaction 

Before taking on the full problem we study the system in 
the absence of Rashba spin-orbit interaction, i.e. for a = 0. 
In Fig. [2] we show the low-energy spectrum and the supercon- 
ducting order parameter when T4c — 5.36 and n = 1. A clear 
superconducting phase transition can be seen between the re- 
gions labeled I U I' and III. Inside the superconducting phase 
I U r there is a sudden jump in the lowest energy levels and 
the appearance of a zero energy state marking the transition 
between I and I'. The energy level jump is correlated with 
an abrupt change in the superconducting order parameter in 




FIG. 2: (Color online) (a) Low-energy spectrum (thin lines), together 
with Zeeman splitting j Vz | (black dotted) and the average of the su- 
perconducting order parameter jA| (thick blue) as function of Vz- 
The two energy levels closest to zero are marked in red. A supercon- 
ducting phase transition can be seen between the F and III regions, 
(b) Real space profile of |A| and (uz). The profile changes con- 
tinuously with Vz, except at the boundary between I and F where 
it changes discontinuously as shown in the figures. Here a — 0, 
Vsc ~ 5.36, n = 1. 



the vortex core, as well as an emergent magnetization of the 
core. In region I the vortex profile is narrow and there is no 
magnetization, while it is much wider in region I' where the 
core is also magnetized. There is thus a competition between 
A and V^ in the core with a local transition taking place at 
the boundary between I and I' . This competition takes place 
throughout the sample and is ultimately responsible for the 
transition into the non-superconducting phase III. However, 
in the vortex core the Zeeman splitting becomes dominant 
earlier because here the superconducting state is suppressed 
by a large rotational component, which gives an additional ki- 
netic energy contribution. Apart from the boundary between 
I and r, the vortex core profile changes smoothly with V^. 
The changes of the vortex core at the transition between I and 
r should be possible to detect experimentally using scanning 
probes. 



B. Including Rashba interaction 

We now continue to study the full Hamiltonian with a ^ 0. 
Figure [3] shows the low-energy spectrum and the supercon- 
ducting order parameter for a — 0.56, Vsc = 5.36, and n = I. 
Two main effects are immediately visible when this result is 
compared to that of Fig. l2] First of all we see a clear asym- 




FIG. 3: (Color online) Similar to Fig. [2] but for spin-orbit coupling 
a — 0.56. The main differences are the appearance of a region 
labeled II and the asymmetry with regards to Vz ■ 



metry in both the energy spectrum and vortex structure as a 
function of Vz- Most notably we only find the I' region for 
positive Vz- We will discuss this asymmetry in depth in Sec- 
tion IIII El The second observation is that the Rashba interac- 
tion has a large effect on the superconducting state, helping 
the superconducting order parameter to survive to higher Zee- 
man splitting. In addition, and for our purposes even more 
important, a finite a creates a smooth transition between the 
region with a large superconducting order parameter to that 
with no superconductivity. This means that a topologically 
non-trivial region II now appears, for which Eq. (Hh is sat- 
isfied. The lowest energy states in region II are Majorana 
fermions, as will be discussed in Section IIIC In Fig. HI we 
show a direct comparison between the superconducting order 
parameter for a — 0.32 and a — 0.48, which explicitly shows 
that a larger a creates a larger topologically non-trivial region 
II. 



In Fig.lSlthe different phases I, I', II and III is visualized in 
{Vz, Vsc) phase diagrams. We once again see the occurrence 
of the r region for positive Vz, but not for negative values. 
The r region occurs as a wedge between I and II for positive 
Vz- In Fig.l6]we plot the low -energy spectrum for the solid 
horizontal cut in Fig. ISljust below where the I' region disap- 
pears and we clearly see that the I' region is not present. On 
the right hand side of Fig. B] we display the {Vz, Vsc) phase 
diagram for a smaller a, showing that the I' region now occurs 
at lower Vz and Vsc- This is to be expected since the Rashba 
interaction helps superconductivity in the presence of a finite 
Vz- Thus, for decreasing Rashba coupling the vortex core gets 
magnetized, i.e. enter region I', at lower values of Vz- 
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FIG. 4: (Color online) Zeeman splitting | Vz \ and average of the su- 
perconducting order parameter j A| for spin-orbit coupling a = 0.32 
(bottom red) and a — 0.48 (top blue). Shaded areas represent the 
topologically non-trivial region II. 
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FIG. 7: (Color online) (a) Probability density for the two lowest en- 
ergy states in Fig. |3] i.e. the Majorana states, (b) Left pair shows 
the same states from a top view, while the three next pairs are the 
probability densities for the next three pairs of states closest to zero 
energy. Note how only the Majorana fermion pair can be separated 
into one core and one edge state. 



C. Majorana fermions 




FIG. 5: (Color online) (14,T4c) phase diagrams for (a) a — 1.04 
and (b) a — 0.8, showing the occurrence of the four regions I, F, II 
and III. The white region has not been mapped. The lines are hand 
drawn approximations drawn on top of data obtained from calcula- 
tions on a 26 X 26 lattice. Solid line marks 14 cut in Fig. ([6|. Dashed 
line marks end point of F wedge for a = 1.04. 
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FIG. 6: (Color online) Similar to Fig. [2] but for spin-orbit coupling 
a = 1.04 (horizontal cut in Fig.Bl. 



So far we have only discussed the topologically trivial and 
non-trivial regions, the superconducting phase transition and 
the vortex core magnetization. While interesting in them- 
selves, our main reason for doing so was to be better equipped 
to identify which features are related to the existence of Majo- 
rana fermions, and which features are due to other effects. We 
now turn our attention to the Majorana fermions which have 
been predicted to occur inside the topologically non-trivial 
phase as zero energy modes in vortex cores.^' ^^''^ Majorana 
fermions, being essentially half of an ordinary fermion, nec- 
essarily comes in pairs, generally located in different vortex 
cores. In a system with only one vortex, the second Majorana 
fermion will instead be located on the edge of the system. In 
Fig. [3] we clearly see two near zero energy modes in the topo- 
logically non-trivial phase II. However, they show slightly os- 
cillatory behavior, such that the mathematical criterion _E = 
for a Majorana fermion in a superconductor is not strictly sat- 
isfied, and there is also several other low-lying states close to 
zero in these regions. 

To establish the Majorana fermion nature of the two low- 
est energy modes in region II we plot the spatial probabil- 
ity density for the i? ~ eigenstates inside region II. If 
plotted directly, these eigenstates appear as two states that 
have finite weight both at the edge and in the core. This is 
an artifact of the ambiguity in basis used to describe these 
(nearly) degenerate eigenstates. In Fig. |7] we therefore plot 
7i = 73 (7+ + 7-) and 72 = ^ (7+ - 7-) rather than 7+ 
and 7_, where 7+ and 7_ are the eigenstates from the numeri- 
cal diagonalization. In the same figure we also plot the proba- 
bility density for the next three pairs of states closest to £■ = 0. 
It is clear that the lowest-energy pair of states (leftmost figure) 
is veiy different from the other states. The two lowest energy 
eigenstates are spatially separable into one core and one edge 
state, clearly displaying their Majorana nature. Distinctly dif- 
ferent, the next lowest eigenstates are non-separable vortex- 
centered states, thus combinable into one single normal elec- 
tronic excitation. These are the equivalent of the Caroli-de 



FIG. 8: (Color online) Same as Fig.JTlbut immediately inside region 
r at the vortex core transition of Fig. [3] The lowest energy pair 
can not be separated into one core and one edge state and is not a 
Majorana pair. 



Gennes-Matricon vortex states found in conventional s-wave 
superconductors . ^* 

Having established the Majorana fermion nature of the 
lowest-energy modes in region II we note that their non-zero 
energies are an artifact of a finite lattice. In fact, it has been 
shown that Majorana fermions that co me close to each other 
experience interference effects.l^SESMOl We therefore draw the 
conclusion that the slight deviations from i? = are due to 
interference between the Majorana fermions in the core and 
on the edge. 



D. Zero mode at the vortex core transition 

Revisiting the vortex core transition that occurs between I 
and r in Fig. [3] we may now ask whether the zero mode that 
appears at this transition is also a Majorana fermion. An im- 
mediate objection to this would be that this is well inside the 
topologically trivial phase, where Majorana fermions not are 
predicted to exist. However, the classification of the topo- 
logical phases are done in a bulk calculation where the or- 



der parameter is constant, while vortex calculations inherently 
implies an order parameter with both varying phase and am- 
plitude. We therefore push this question a little further. In 
Fig. (jHll we plot the probability density for these states. It 
is clear that the first pair of states do not separate into one 
core and one edge state as they did in Fig. |7] and we can con- 
clude that this zero mode is a normal fermionic state, and thus 
very different from the Majorana fermion that exists inside the 
topologically non-trivial phase II. This zero-energy state thus 
presents an experimental challenge when searching for Majo- 
rana fermions. 



E. Spectral asymmetries 

Finally, we address the asymmetry that is seen in the en- 
ergy spectrum when a ^ 0. It is clear from our numerical re- 
sults that there is an asymmetry in the energy spectrum as well 
as vortex core order parameter and magnetization as function 
of Vz- In Fig. l9]the result of flipping the vortex direction 
n —>■ —n and the sign of the spin-orbit interaction a —> —a is 
shown for \a\ — 0.48. The energy spectrum is independent of 
the sign of the spin-orbit coupling, but the dependence of Vz 
is flipped for reversed vortex rotation. The energy spectrum is 
thus asymmetric with respect to Vz — > —Vz and n — 7> —n, but 
the results shows a symmetry with respect to the simultaneous 
parameter transformation V^, n — >■ —Vz, —n. 

In order to gain more understanding of the physical origin 
of these symmetries and asymmetries we consider the analyt- 
ical continuum equivalent of Eq. ([T}. The Hamiltonian, in the 
gauge where A has a constant phase, can then be written as^ 



H = 



D1 + D2 
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Noting that the numerical calculations indicate that the sys- 



tem is symmetric for a = 0, we start our analysis in this 
limit. The S terms then disappears and it is easy to see that 
Vz —> —Vz implies D3 —> —D3, while all other terms are 
unchanged. The effect on the Hamiltonian of this physical 
transformation, which corresponds to a reversal of the Zee- 
man field direction, can almost be undone by the similarity 
transformation sHs^^, where 



(7) 



This transformation corresponds to a relabeling of up and 
down spins, or equivalently, a spin flip. The only place sHs^^ 
differs from H is a minus sign in front of the A term, which 
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FIG. 9: (Color online) Similar to Fig. plbut for spin-orbit coupling \a\ — 0.48 and different signs on a and vortex rotation n. The energy 
spectrum is asymmetric with respect to V^ — >■ — 14 and n — >■ —n, but symmetric with respect to a — >■ —a and T4 , n — >■ — 14 , —n. 



does not affect the energy spectrum. Physically, the Zeeman 
term breaks the symmetry along the z-axis, and when it is 
reversed it suffices for any eigenstate in the original system 
to reverse its spin to become an eigenstate with the same en- 
ergy in the new system where the Zeeman field is reversed. 
Continuing our analysis at a == we also look at the effect 
of n — > —n. This transformation, which corresponds to a 
reversal of the vortex rotation direction, gives D2 -^ — -D2- 
This physical transformation can be undone by the coordinate 
transformation 6 ^> —6. This means that any eigenstate in the 
system with vortex rotation direction n can be brought into 
an eigenstate with the same energy for a system with vortex 
rotation —n by simply reversing the rotation direction for the 
state. 

The situation becomes more complicated for finite a. The 
transformation V^ -^ —Vz can no longer be undone by the 
similarity transformation sHs^^, because the off-diagonal S 
elements are interchanged during this process. Similarly, the 
physical transformation n — > —n is no longer counteracted 
by the coordinate transformation 9 — )■ —9, because, again, 
the S terms are modified in a non-trivial way in this process. 



However, if we reverse both the magnetic field and the vortex 
direction, i.e. Vz,n —> —Vz,—n, then the coordinate trans- 
formation 9 —i' —9, followed by the similarity transformation 
sHs^^ brings back the original H apart from an additional 
minus sign on all off-diagonal terms, which does not influ- 
ence the energy spectrum. 

Mathematically, the spectral asymmetries in Vz and n are 
directly related to the introduction of the Rashba spin-orbit 
dependent S terms in Eq. dSll. To understand this we note that 
when only the Rashba interaction is present, such that the only 
term that is non-zero is S^., then the spin lies in-plane and is 
locked perpendicular to the direction of motion. This is evi- 
dent from the Rashba spin-orbit interaction in Cartesian coor- 
dinates: a(k X (t)z. It is also clear that the energy depends 
on which of the two directions perpendicular to the direction 
of motion the spin locks onto, and that the energy is preserved 
under a simultaneous reversal of the directions of motion and 
spin. Coming back to the full problem, the Zeeman term will 
cause an out-of-plane tilting of the spin. Taking also into ac- 
count the kinetic and superconducting terms, the exact form 
of the eigenstates can be expected to be quite complicated. 



However, for our purpose it is not important to know the exact 
eigenstates of the Hamihonian, rather it suffices to draw the 
following conclusions about how three of the terms influence 
the eigenstates: 

• D3: Each eigenstate gets an energy contribution due to 
the Vz term, which depends on the spin being tilted up- 
or downwards. 

• Sj. (Sj. ): Each eigenstate gets an energy contribution 
from the Rashba spin-orbit interaction, which depends 
on the relative orientation of the direction of motion and 
the in-plane spin component. 

(2) 

• D2 (Sj. ): Each eigenstate gets an energy contribution 

from the superconducting vortex term, which depends 
on the particle motion being with or against the rotation 
of the vortex. 

(2) 
Above Sj. appears in parentheses because it contributes to a 

mixture of the two latter effects. With regards to the spin-flip 

operation, the ±-index transforms in the same way as for the 

However, under the n -^ —n and — > —9 trans- 



S-^^ term. 



:(2 



formations, S±. is multiplied by a negative sign relative to the 

Sj. term. With this we are finally able to explain the results 
in Fig. I9] We begin with the observation that for Vz -^ —Vz 
the off -diagonal S terms obstructs attempts to get back eigen- 
states with the same energy through a spin flip. The reason 
is that the Rashba term does not in general contribute the 
same amount of energy to the eigenstates after this transfor- 
mation, unless the rotation direction for the states are reversed 
at the same time, which requires the coordinate transforma- 
tion 6 —i' —6. However, the energy contribution from having a 
particle rotating either clockwise or counter-clockwise around 
the vortex is not the same unless the vortex rotation direction 
is also changed. Therefore the energy spectrum is preserved 
only under a simultaneous reversal of Vz and n. A schematic 
depiction of this analysis is given in Fig. 



10 We see that the 



asymmetry in the energy spectrum is the combined effect of 
Zeeman, Rashba, and vortex rotation. The n -^ —n asymme- 
try can be worked out following a completely analogues pro- 
cedure and can be visualized by reading the steps in Fig.[TO]in 
reversed order 

Finally, we note that the symmetry observed under the re- 
versal of a is very similar to that expected from a global phase 
transformation in A. This can be seen by concluding that a 



-5« 



-4^). 



—a has the effect of transforming 3^,3^. 

Such a transformation leaves det{H — IE) invariant and the 

energy spectrum is therefore guaranteed to be preserved. 



IV. CONCLUSIONS 

We have studied a ferromagnetic 2D Rashba spin- 
orbit coupled semiconductor with proximity-induced spin- 



singlet .5-wave superconductivity by using the self-consistent 
Bogoliubov-de Gennes method on a square lattice. The calcu- 
lations confirm the existence of Majorana fermions in vortex 
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FIG. 10: (Color online) Schematic figure depicting the effect of ap- 
plying parameter and coordinate transformations to the Hamiltonian 
in Eq. Jsl. Only the relative sign between the S terms is important for 
the energy spectrum: det{H — IE) — 0, and the first and last stage 
is therefore energetically equivalent. The three middle stages are al- 
ways energetically equivalent as they only differ by a (non-physical) 
coordinate transformation. 



cores in the topologically non-trivial phase. Strong Rashba 
spin-orbit interaction makes the non-trivial phase larger as it 
counteracts the pair-breaking effects of the Zeeman spin split- 
ting on the superconducting state. We have also established 
the existence of a second superconducting phase transition 
within the topologically trivial phase, characteiized by a fi- 
nite vortex core magnetization and a wider vortex profile. A 
normal fermionic zero-energy mode appears at the transition 
between these two superconducting states. This demonstrates 
the existence of non-Majorana zero-energy excitations in the 
vortex core. We have also found a pronounced asymmetry 
in the vortex core energy spectrum with respect to the Zee- 
man splitting and the vortex rotation direction. This asymme- 
try is a consequence of the interaction between the, partially 
competing, Zeeman, Rashba, and vortex rotation terms in the 
Hamiltonian, resulting in energy states with a delicate depen- 
dence on the spin and orbital coordinates. In aggregate, these 
results show on a complex behavior of the vortex core energy 
spectrum in Rashba spin-orbit coupled superconductors. 
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